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Abstract. Using the algebraic dynamical method, the dynamics of entanglement in 
an atom-field bipartite system in a mixed state is investigated. The atomic center-of- 
mass motion and the field-mode structure are also included in this system. We find 
that the larger values of the detuning and the average photon number, the smaller 
fluctuation of the entanglement, but the period for the evolution of the entanglement 
doesn't increase accordingly; It is also found that the fluctuation of the entanglement 
varies slightly with the atomic motion and oscillates fast with the value of the field- 
mode structure and the transition photon number increasing. Moreover, a damping 
evolution of the entanglement appears when considering detuning and the atomic 
motion simultaneously. 
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1. Introduction 

Entanglement has been widely investigated in quantum information processing [U 
|2j |3l HI [5] and plays important roles in many potential applications, such as 
quantum communication, quantum teleportation, quantum cryptography, entanglement 
swapping, dense coding, and quantum computing, etc [6j [7J EJ [9j [10J El [12] • Moreover, 
with the recent rapid developments in quantum information, more precise quantifications 
for entanglement are required to produce, manipulate and detect entanglement states. 
Most of earlier studies pay much attention to the bipartite systems that start in 
pure states, which are presently well understood [13]. While for mixed states, the 
entanglement is difficult to be measured, because it is not easy to define an analog of 
the Schmidt decomposition. However, it is more significant for studying the systems 
in mixed states since there doesn't exist a system that could be decoupled perfectly 
from environmental influences for realistic states observed in experiments. Several 
entanglement measures for mixed states have been proposed, such as the entanglement 
of formation Ep(p), the entanglement cost Ec{p) and the distillable entanglement 
Ed(p) [El US], but these are not effective computational means. To overcome these 
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shortcomings, Peres and Horodecki have proposed a new criterion for separability 
[T6| [T7] . Based on this criterion, negativity is introduced, which is a good quantity 
for measuring the entanglement of mixed states [EJ EE] and is equal to absolute value 
of the sum of negative eigenvalues of the partial transpose of density operator with a 
bipartite mixed state. 

There are considerable researches on the dynamics of entanglement for the Jaynes- 
Cummings model (JCM) [20J with a mixed state. The entanglement for the model 
consisting of a two-level atom interacting with a single mode field through a two- 
photon process was studied [21]. In [22], the authors devoted to studying the effects 
of phase damping on the entanglement evolution, where the two-level atom initially 
in a mixed state interacts with a multi-quanta single mode quantized cavity field 
including cavity loss. Akhtarshenas and Farsi |23j studied the entanglement of the 
JCM with the atom initially in mixed state and the field prepared in Fock states. 
Very recently, the entanglement for the JCM with a two-level atom in a lossy cavity 
was manipulated by a classical driving field [21], which was a novel scheme to control 
the behavior of the entanglement. Xu and co-workers [25J studied the entanglement 
dynamics of a three-level atom in H configuration, and they found that the maximal 
value of entanglement decreases with the temperature of the thermal field and can 
be controlled by the detuning. Zhang and Xu [26] considered the entanglement of 
JCM interacting with two mode quantized fields for two different initial situations. 
Yonac [27] et al studied the entanglement dynamics of two separate Jaynes-Cummings 
Hamiltonians which consist of two identical two-level atoms inside two separate cavities, 
respectively. They investigated how the two sets of pure initial states affected the 
entanglement dynamics with resonant situations. The literature [28j has shown that the 
pure entangled states are very difficult to obtain due to the interaction of uncontrollable 
environment and thus the atomic initial state prepared in mixed states is more practical 
and general. In this paper, we will focus on the problem of a moving two-level atom 
coupled with a quantized field in a multi-photon transition process. The effects of the 
detuning, the transition photon number, atomic motion, the field-mode structure and 
the temperatures on the entanglement are all discussed. In fact, the model including all 
these factors is extremely difficult to deal with using the conventional methods. Here 
we adopt the algebraic dynamical approach [291 ED] to study the nonlinear system, and 
the key idea is introducing a canonical transformation that transforms the Hamiltonian 
into a liner function in terms of a set of Lie algebraic generators. The nonlinear system 
is thus integrable and solvable, then the density operator with time evolution can be 
derived easily. 

The paper is organized as follows: In Section 2, we present the basic model, and 
derive the time evolution expression of the density matrix using the algebraic dynamical 
approach. In Section 3, we study the dynamics evolution of the entanglement between 
the two subsystems for different situations with the atom initially in a pure state and a 
mixed state. Finally, a brief conclusion is given in Section 4. 
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2. The multi-photon JCM 

We investigate a moving two-level atom interacting with a single mode quantized 
radiation field in a multi-photon transition process. Within the rotating wave 
approximation, the Hamiltonian of the bipartite system can be described as 

H = uj S z + uja)a + g\f(z)}\a) l S^ + a l S + ) {h=l), (1) 

where Uq is the atomic transition frequency, u is the field frequency (a) denote the 
creation (annihilation) operators of the radiation field and S + (S-) represent the atomic 
raising (lowering) operators. S z is the atomic inversion operator, g is the atom-field 
coupling constant. f(z) is the shape function of the cavity field mode [31] and the 
atomic motion can be incorporated as 

/(*) f(vt), (2) 
where v is the atomic motion velocity. A specific mode TEM mnp is defined: 

/W = sm(^), (3) 

with p represents the number of half-wavelengths of the field mode inside a cavity of the 
length L. When the atomic motion velocity is considered as v = gL/ir, the evolution 
function 9(t) can be obtain 

9{t) = [\f(vt')] l dt>. (4) 



We restrict the motion of the atom along the z-axis. There exists one conservative 
quantity for the system as follow: 

N = a?a + l(S z + ~), (5) 

which commutes with Hamiltonian (1). To linearize the Hamiltonian by the algebraic 
dynamical approach, a canonical transformation operator is introduced as 

U a = ev[-£ 75 {a l S + -a*S-)], (6) 

where a = - arctan[(^/A 2 /4 + g 2 F - A/2)/gF 1/2 }, A = u - kj, g = g9(t)/t, and 
F = N\/(N — Based on the algebraic dynamical method, SU(2) algebra generators 
{J , J+, J-} are introduced, with J = S z , J + = F~ l / 2 a l S + , J_ = F^^a^S^, which 
satisfy the following commutation relations: 

[Jo, J+) = J+, [Jo, J-] = -J-, [J+, J-] = 2J - (7) 
Then the canonical transformation operator is reduced as 

U g = exp[a( J+-J-)]. (8) 

After the canonical transformation, the dressed Hamiltonian can be written in terms of 
the SU(2) algebra generators 

H' = U-\ [ Hdt)U g = E(N) + A J , (9) 
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A = y/A 2 + Ag 2 F, E(N) =u(N -1/2). 
Then the time evolution operator can be expressed by 

u(t) = e -^ Hdt 

= U g exp[-i(E(N) + \J )]U g 1 

= e- iE{N)t [cos y - 2tJ sin y cos 2a + i(J + + J J) sin y sin 2a}. (10) 
Assuming the cavity field is initially in the single-mode thermal state 

oo n 
n=0 ^ ' 

with m = l/[exp(o;/T) — 1], is the mean photon number of the cavity field for thermal 
equilibrium at a certain temperature T, and the atom is initially in a statistical mixed 
state 

Pa (0) = C e \e)(e\ + C g \g)(g\, (12) 

with C e +C g =l, C g = (C g = l) means the atom is initially prepared in the excited 
state (the ground state) and < C g < 1 represents the initial state of the atom is 
in the mixed state. Since there is no system that could be decoupled perfectly from 
environmental influences and thus mixing is unavoidable, mixed states are more realistic 
and meaningful than pure states in experimental situation [32]. Accordingly, the initial 
state of the total system can be given as 

Pa/(0) = p/(0) g> p a (0) = PnC g \n, g)(n,g\ + J2 PnC e \n, e)(n, e\, (13) 

n n 

Using Eqs.(10) and (13), the density operator with the time evolution can be written 
by 

Paf (t) = U(t) Paf (0)U+(t) = Y,PnC 9Pg (t) + J2PnCe P e(t), (14) 



where 



p g (t) = [cos 2 (^) + sin 2 (^) cos 2 (2a n )]\n,g)(n,g\ 

— i sin(— — ) sin(2a n )[cos + i sin cos(2a n )]|ri, g)(n — I, e\ 

+ i sin(^-) sin(2a n )[cos — i sin cos(2a n )]|n — /, e)(n,g\ 

2 2 2 

+ sin 2 (-^-) sin 2 (2a n )|n — I, e)(n — /, e|, (15) 
p e (t) = [cos 2 (^y~) + sin 2 (^y^) cos 2 (2a n+ i)]\n,e)(n,e\ 

■ . /A n+ ;t, . , u X n +lt . . X n+ it , ,,. . . . 

- % sm(— — ) sm(2a n+i ) [cos — - i sin — — cos(2a„ +i )J \n, e){n + l,g\ 

+ i sm(^Y~) sin(2a n+i ) [cos h^L + i s j n JH^l. CO s(2a n+z )] |n + I, g) (n, e| 
+ sin 2 (^±^) sin 2 (2a n+i )|n + l,g)(n + l,g\, (16) 
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with A n = v/A 2 + Ag 2 F(n) , a n = — arctan[(A/A 2 /4 + g 2 F(n) - A/2)/(gF 1 / 2 (n))], 
F(n) =n\/(n-l)\. 

3. Entanglement evolution for the bipartite system 

In this section negativity is employed to investigate the entanglement for the coupled 
atom-field system. The system is described by the density matrix in Eq.(14), and the 
entanglement degree is quantified by [18] 



Af(p) = (IT) 

where p T is the partial transpose of p. ||p T || denotes the trace norm of p T which is 
equal to the sum of the absolute values of all eigenvalues of p T . M{p) is equal to the 
absolute value of the sum of the negative eigenvalues of p T , which has been proposed 
as a quantitative entanglement measure. It is a good measure of entanglement because 
negativity is an entanglement monotone, does not increase under local operations and 
classical communication (LOCC) and can also be easily treated mathematically. From 
the final state in Eq.(14) and using Eq.(17), the amount of entanglement can be obtained 
from the following result 

l—l oo 



\P T 



2 

n=0 n=0 



+ \Hn + £n+l - a/ (Pn ~ £,ri+l) 2 + 



iXr. 



where 



Hn = C 9 P n [cos 2 (^-) + sin 2 (^-) cos 2 (2a n )] + C e P n _ { sin 2 (^-) sin 2 (2a n ), 

^ n = C e P n [cos 2 (^-) + sin 2 (^-) cos 2 (2a n+ i)] + C g P n+i sin 2 (^±^) sin 2 (2a n+i ) 



n 



= i(C g P n+ i - C e P n ) sin(^— ) sin(2a n+i )[cos i sin cos(2a n+ |)],(19) 

Based on Eqs.(17) and (18), we present some interesting numerical results for 
different parameters to demonstrate the effects on the entanglement degree. In Fig. 
1, the negativity is shown as functions of C g and the scaled time gt/n. When the 
probability of the atom initially in the ground state is lager, the entanglement is smaller, 
and the entanglement is close to zero as the value of C g exceeds 0.5. Therefore, we 
choose C g = 0.2 to demonstrate the effects of the parameters on the entanglement. 
Fig. 2 displays the relation between the entanglement and the atomic motion and 
the field- mode structure. Here, the atomic motion is considered as v = gL/ir and 
is neglected by performing 9{t) — > t. It is found that the atomic motion leads to the 
regular evolution period of the negativity by comparing Fig. 2(a) and Fig. 2(b), and the 
increasing of p shortens the period of the entanglement and results in a slight decreasing 
of the amplitude in Figs. 2(b) and (c) (In fact, when the parameter varies in a small 
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region, the amplitude of the entanglement doesn't change a lot). Interestingly, when 
the atom motion is considered and the evolution time is in certain times gt = 2nn/p 
(n = 0, 1, 2 • • •), the negativity goes to zero which is called "entanglement sudden death" 
(EDS) [331 EH [35] and in the other Figs, the similar phenomenon also occurs. 

Now we shall turn our attention to the effects of the detuning. What is worth 
noting that most of previous researches only consider the case of the resonance to 
avoid the complexity. Without loss of generality, considering the detuning is necessary 
and meaningful. In Fig. 3, it clearly shows that when the detuning increases the 
entanglement decreases and oscillates fast [36], and it is interesting to note that when 
the detuning and the atomic motion are all considered, a damping evolution of the 
entanglement arises. The former is attributed to the fact that on having large dutunings, 
the atomic system is weakly coupled to the radiation field, and the latter may be referred 
to a quantum decoherence [371 E3 [39] when the interaction between the internal and the 
external atomic dynamics is considered (see Fig.3(b)). Also, it is shown that the atomic 
motion can lead to a damping in the correlation functions of non classical behaviors 
between the field and internal atomic variables, which induces the separability of the 
atom and the field in [38]. Fig. 4 depicts how the average photon number m influences the 
entanglement evolution. One can see that with increasing the average photon number 
the negativity decreases, but the evolution periodicity of the entanglement doesn't 
change. These results show that the maximum value of the entanglement depends on 
the mean photon number of the field (or the thermal equilibrium temperature T), but 
the period doesn't, which agrees well with the results in [10] • As a matter of fact, for the 
thermal field, when the average photon number increases, the weight mixture of number 
states denoted by P n becomes small, so the entanglement feature is washed out greatly. 
In addition, different photon transition processes can also affect the entanglement, as 
is shown in Fig.5. Intensity of the entanglement increases as the photon number I 
increases, meanwhile the oscillation of the negativity becomes faster and faster [4"T] . 

Experimentally, the realization of the JCM in this paper is not difficult. We could 
inject an highly excited Rydberg atom along the axis of the cavity which moves fast to 
eliminate no reflection effect from the cavity, then we could assume that it has a constant 
velocity along the cavity [42]. The transition between the upper and lower levels of 
the atom may involve I photons by adjusting the energy separation between the level 
close to the energy of I quanta of the electromagnetic field. Photons being the fastest 
information carriers have only very weak coupling to the environment, and are thus the 
most successful physical system for the observation of multi-partite entanglement and 
proof-of-principle demonstrations of a diversity of quantum information applications so 
far. Multi-photon entanglement has become a hot topic recently and we believe the new 
approach in this paper will promote the relate development. 
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4. Conclusion 

We have used negativity to study the entanglement for the system of a moving two-level 
atom interacting with a quantized single mode field in the mixed state with a multi- 
photon process. With the algebraic dynamical method, its solution is obtained. In 
our work, the atomic motion and the field-mode structure, the detuning, the transition 
photon number and the average photon number of the field are all considered. Our 
studies show that increasing of the detuning or the average photon number can lead 
to decreasing of the negativity, the entanglement oscillates fast for lager detunings and 
EDS appears with certain times. Atomic motion and the field-mode structure also give 
rise to many interesting effects on the entanglement. An increase of the parameter p 
can shorten of the evolution periodicity of the entanglement. Importantly, a damping 
evolution of the entanglement happens when the detuning and the atomic motion are 
considered together. Furthermore, we demonstrate the fact that a lager photon number 
/ makes it quick to the oscillation of the negativity and bring its amplitude up slightly. 
Our researches might shed light on understanding the entanglement of the JCM and 
provide a simple means to study dynamics of the entanglement. 
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Figure 1. The evolution of the negativity as functions of C g and the scaled time gt/n. 
The atomic motion is considered, other parameters: m — 1, A = 0, p = 1, 1 = 2. 
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(c) 

Figure 2. The evolution of the negativity as a function of the scaled time gt/w for 
different parameters p. Other parameters: m = 1, A = 0, I = 1, C g = 0.2, (a) the 
atomic motion is neglected, (b) p — 1, (c) p = 3. 
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(a) (b) 

Figure 3. The evolution of the negativity as functions of the scaled time gt/ir and A. 
Other parameters: m = 1, I = 2, C g = 0.2, (a) the atomic motion is neglected, (b) the 
atomic motion is considered (p — 1). 




Figure 4. The evolution of the negativity as functions of the scaled time gt/ir and 
the average photon number m. The atomic motion is considered, other parameters: 
A = 0, p= 1, I = 2, C g = 0.2. 
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Figure 5. The evolution of the negativity as a function of the scaled time gt/w for 
different parameters I. The atomic motion is considered, other parameters: m = 1, 
A = 0, p = 1, C g = 0.2, (a) I = 1, (b) I = 3, (c) I = 8. 



